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ABSTRACT 

Cyclical period changes are a fairly common phenomenon in close binary systems and 
are usually expla ined as due to either the magnetic activity of one or both components (e.g., 
lApplegatel 1 1 992l) or to the light-travel time effect(LTTE) of a third body. W e searched th e 
orbital period changes in 1 82 EA-type (including the 101 Algol systems used by lHall I (ll989l) X 
43 EB-type and 53 EW-type binaries with known both the mass ratio and the spectral type of 
their secondary components. We reproduced and improved the same diagram as Hall's ( 1989) 
accord ing to the new collected data. Our plots do not support the conclusion derived by iHall I 
(Il989l) that all cases of cyclical period changes are restricted to binaries having the secondary 
component with spectral types later than F5. The presence of period changes also among stars 
with secondary component of early type indicates that the magnetic activity is one cause, but 
not the only one, for the period variation. It is discovered that cyclic period changes, likely due 
to the presence of a third body are more frequent in EW-type binaries among close binaries. 
Therefore, the most plausible explanation of the cyclical period changes is the LTTE via the 
presence of a third body. By using the century-long historical record of the times of light 
minimum, we analyzed the cyclical period change in the Algol binary WW Dra. It is found 
that the orbital period of the binary shows a ~ 112. 2yr cyclic variation with an amplitude 
of ~ 0.1977days. The cyclic oscillation can be attributed to the LTTE via a third body with 
a mass no less than 6.43M0. However, no spectral lines of the third body were discovered 
indicating that it may be a candidate black hole. The third body is orbiting the binary at 
a distance shorter than 14.4 AU and it may play an important role in the evolution of this 
system. 

Key words: binaries : close - binaries : eclipsing - stars: individual: WW Dra- stars: late- 
type. 



1 INTRODUCTION 

Orbital period changes of stellar eclipsing binary systems can be 
investigated by analyzing the (O - C) diagram showing the dif- 
ference between the observed epochs of light minimum and those 
computed with a given ephemeris. A periodic pattern in an O - C 
curve is a fairly common phenomenon in Algols, W Ursae Ma- 
joris bina ries, and RS Canum Venaticorum and the cataclysmic 
variables llHalll 1 19891 : iHall & KreineJ 1 19801 : iHobart et al. 1 1 1994 
IWamerlll988r) . Similar patterns of the O - C diagram for several 
classes of close bina ries suggest a common underlying mechanism 
jZavala et al. 1 120021. such as mass loss, apsidal motion, magnetic 
activity, and presence of a third body. iZavala et al. 1 12002) thought 
that apsidal motion and mass loss are unlikely mechanisms. There- 
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fore, at present, cyclical period changes can usually be explained 
as du e to either the m agnetic activity of one or both components 
(e.g.. lApDlegatelll992h or to the light-travel time efFect(LTTE) via 
the presence of a third body. 

The hypothesis that cyclical period changes are caused by 
the presence of a third body has been discussed by several in- 
vestigators ( Frieboes-Conde & Herczeg 1973; Chambliss 1992J; 
iBorkovits & Hegedils Ill99q) . In this hypothesis, the motion of the 
binary around the center of mass of the system causes a peri- 
odic change in the observed period due to a light-travel time ef- 
fect(LTTE), thereby creating a periodic pattern in O - C curve. 
Afterwards, Hall found a striking correlation between the spectral 
type of the low-mass secondary com ponent and the presence of 
a cyclical period change in his study l lHalllll98^) on 101 Algols. 
From his plot, he noted that all cases of cyclical period changes 
are restricted to systems with spectral types of the secondaries later 
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than F 5. Based on this result. lApplegate I ( Il992n and iLanza et al. I 
( 119980 developed a theory to explain the periodic pattern in O - C 
curves of these systems. In this theory, a certain amount of an- 
gular momentum is periodically exchanged between the inner and 
the outer parts of the convection zone, and therefore the rotational 
oblateness of the star and hence the orbital period changes while the 
system's component goes through its acti vity cycles. However, the 
period changes of those Algols used by Hall ( 1989) were mainly 
derived from visual and photographic observations. This theory 
was frequently used to interpret the orbital period modulation 
of close binaries containing at least one cool component (e.g., 
[Hall 1991; Oian et al. 1999; Oian et al. 2000). As new and more 
accurate observational material has accumulated since then, in the 
present work, we will reproduce and improve the same diagram as 
Hall's (1989) of EA, EB and EW-type binaries based on the new 
collected data. We will check Hall's plot and discuss the cause of 
cyclical period changes. Meanwhile, we will analyze the cyclical 
period change of the RS CVn-type binary WW Dra derived from 
the century-long historical record of the times of light minimum 
and discuss its plausible cause. 



2 THE MOST PLAUSIBLE EXPLANATION OF THE 
CYCLICAL PERIOD CHANGES 

As discussed above, at present, the m agnetic activity of one or 
both components (e.g.. lApDlegatell992i) or the light-travel time ef- 
fect(LTTE) via the presence of a third body are usually invoked 
to exp lain the cyclical period changes of close binaries. Hall I 
( Il989f ) searched the orbital period changes of 101 Algol-type bina- 



ries in iGiuricin et eJ. 1 1 19830. The samples of our study are made 
by the stars listed in iKreiner et al. I ( 1200 ih , the 101 Algol sys- 



tems in|Giuricin et al. I ( Il983l) . and the Algol-type binaries listed in 
libanoglu et al. I J200a ). As selection criterion we considered stars 
that either show cyclical period changes or have secondary compo- 
nent of late spectral ty pe. Fi nally, 182 EA-type (including the 101 
Algol systems used bv lHalll |l989)), 43 EB-type, and 53 EW-type 
binaries were collected for this study. 

In this paper, the data of EA, EB, and EW-type binaries are 
presented in Tables 1-3, respectively. In Table 1, Column (I) and 
(7) give the systems we selected; (2) and (8) the secondary com- 
ponent's spectral type; (3) and (9) the mass ratio; (4) and (10) the 
form of the period change; (5) and (II) the geometrical structure 
of the binary : semidetached binaries(SD) or detached binaries(D), 
and (6) and (12) the reference for the O - C information. 

In Table 2 and 3, Column (1) gives the systems we selected; 
(2) the secondary component's spectral type; (3) the mass ratio; (4) 
the form of the period change and (5) the reference for the O- C in- 
formation. The secondary component's spectral type (Sp,) and the 
mass ratio {q) are up to date values taken from one of these refer- 
ences: iKreiner et al. I f200lh . lGiuricin et al. l ( ll983h . libanoglu et al. I 
(;2006j, the reference given in corresponding table, and VizieR 
database[J. Therefore, some of the secondary component's spectral 
types in Table 1 differ from those in lOiuricin et al. I ( 11983' ). More- 
over, in the process of investigation, we reclassifi ed several sys- 
tems a s EB type or detached binaries with respect to lCiuricin et al. I 
( 119831) according to the mentioned more recent bibliography. 

The plots of mass ratio(^) vs. secondary component's spectral 
type (SP2) for EA, EB, and EW-type binaries are displayed in Figs. 
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1-3, respectively. The form of the period change follows the con- 
vention adopted bv lHall I jl989l) . A horizontal line ( - ) indicates no 
period change, a forward slash ( / ) indicates a period increase only, 
a back slash ( \ ) indicates a period decrease only, a cross ( x ) indi- 
cates both increase and decrease of the period, and a filled circle ( • 
) is used for systems for which we have inadequate data for judge- 
ment. In Fig. I, the magenta symbols are used for the semidetached 
Algol-type binaries and the black ones are for detached Algol-type 
binaries. It is clear fr om Fig. I th at our plots do not support the 
conclusion derived bv lHalll d 19891) that all cases of cyclical period 
changes are restricted to binaries with secondary component with 
spectral type later than F5. There are cases among both semide- 
tached and detached Algols in which the spectral type of the sec- 
ondar y compon ent is earlier tha n F5, such as RW Cap (Erdem et al. I 
|2007|) . TX Her JAk et al. |2004|) . and it is expected that the number 
of these systems will grow rapidly as more new observational data 
will be derived. The presence of period changes also among sys- 
tems with low-mass component of early-type stars rules the mag- 
netic activity out as unique cause for the period variation. Among 
binaries with late-type component the orbital period variation can 
be due to either magnetic activity or LTTE. Whereas, among bina- 
ries with early-type components the LTTE is the more likely cause. 
Moreover, the validi ty of the Applegat e mechanism has recently 
come into question ( Lanza 112003 , 120061) . Lanza suggested that the 
Applegate mechanism should be rejected because it can not explain 
the orbital pe riod modulations of classical RS CVn close binaries 
dLanzaluOOSh . Afterwards he also found that the mechanism is in- 
adequate to explain the cyclic al period cha nges of all close bina- 
ries with a late-type secondary ( lLanzal2006l) . Again, the Applegate 
mechanism predicted that there is a connection between the lumi- 
nosity variation and the variation of period. However, to date, no re- 
liable connections were found in the literature. Therefore, the most 
plausible explanation of the cyclical period changes is the LTTE via 
the presence of a third body. We found that 48.9 % of EA, 44.2 % 
of EB and 64.2 % of EW-type binaries have cyclical orbital period 
variation. If we assume that such variations are related to the pres- 
ence of a third body through the LTTE, then we find that EW stars 
have the highest probability to belong to mu ltiple systems. Thes e 
results are in agreement with the findings of IChambliss I dI992bl) . 
The detailed statistical numbers of cyclical period changes in close 
binary systems are displayed in Table 4. 

In the following sections, we present our investigation on the 
cyclical period change in the RS CVn-type binary WW Dra and 
discuss about its causes as the presence of a black hole companion. 



3 NEW CCD PHOTOMETRIC OBSERVATIONS FOR 
WW DRA 

WW Dra ( = HD 150708 = HIP 8I5I9 = BD +60°I691, 
Vmax - 8 . 3 mag ) was discovered to be an eclipsing binary by 
iHarwoodl JI9I6|) . It i s a RS CV n-type eclipsing binary with 
G2+K0 spectral type J.Tov Ill94lh . Studies based on photographic 
and photoelectric observations were carried out by iPlaut I (Il94(]|). 
iMezzetti et^ ( ll979l) . lMardirossian etlTI jl980l) and lTunca et al. I 
1 19811) . Some of them also calculated the orbital and physical ele- 
ments of WW Dra, and the binary was confirmed to be a detached 
system composed of two sub-giant stars. The pe riod variation of 
this binary was studied bv lAlbavrak et al. I d 19991) who derived the 
parameters of the light time orbit. However, as more new observa- 
tional data have been derived since then, we will display different 
results of orbital period change of WW Dra. 



Table 1. Period changes of semidetached and detached Algol-type binaries. 
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Table 2. Period changes of EB-type binaries. 



Table 3. Period clianges of EW-type binaries. 
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In order to analyze the period variations of WW Dra and in- 
vestigate the physical properties of the third body, CCD observa- 
tions were acquired on May 31, 2007 with the PI 1024 TKB CCD 
photometric system attached to the 1.0-m reflector at the Yunnan 
Observatory. The V filter, close to the standard Johnson UBV sys- 
tem, was used. The effective field of view is about 6'. 5 x 6'. 5 at the 
Cassegrain focus and the size of each pixel is 0".38. The integration 
time is 60 s for each image. The coordinates of the nearby compar- 
ison star are RA=16:39:03.91, DEC=-h60:42:02.6 (J2000.0). The 
PHOT task of IRAF, which measures the aperture magnitude for a 
list of stars, was used to reduce the observed images. By using our 
photometric data, we provided the most recent determination of 
time of light minimum - HJD 2454252.3029(±0.0018). 
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Table 4. Statistical numbers of cyclical period changes in close binaiy sys- 
tems. 





EA-type 


EB-type 


EW-type 


(l)TotalNo. 


106 (SD) + 76 (D) 


43 


53 


(2)No. of X 


66 (SD) + 22 (D) 


19 


34 


(3)Binanes with Sp2 earlier than F5 in (2) 


6 (SD) + 7 (D) 


15 


6 


(4)ratio+ 


48.9% 


44.2% 


64,2% 


(5)ratio* 


14.8% 


78.9% 


17,6% 



4 ORBITAL PERIOD VARIATION OF WW DRA 

To investigate the physical properties of the third body in WW Dra, 
we search for cyclic orbital period changes. A total of 92 times of 
light minimum from the literature have been collected and com- 
piled in the present paper Most of times of light mi nimum were re- 
trieved from compilation of Hall & Kreineg (198Cf). Times of light 
minimum are listed in the first and eighth column o f Table 5. In our 
analys is, the value HJD244322 1.544 obtained bv ISudding et al. I 
( 119771) was not used because its (O - C) value shows large scat- 
ter when compared with the general trend formed by the other data 
points. In the second and ninth column the number (E) of orbits 



Note. ^ Ilie ratio of binaries sliow eyeliea] period clianges to tlie total number. 

* the ratio of binaries not only show cyclical period changes but also have Sp2 earlier than F5 to biniu 

period changes. 



elapsed from the initial epoch of primary minimum is listed. The 
{0 — C)i values of all times o f light minimum were computed with 
the linear ephemeris given bv lKreiner et al. I ( 1200 II) . 



Min.I = HJD 2427983.3236 + 4''.6296328£, 



(1) 



They are listed in the third and tenth column of Table 5 and 
plotted vs. number of orbital periods in the upper panel of Fig. 
4, where open circles refer to photographic or visual observations. 
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Figure 1. Plot of mass ratio((;) vs. secondary component's spectral type 
(SP2) for EA-type binaries listed in Table 1. A horizontal line ( - ) indicates 
no period change, a forward slash ( / ) indicates a period increase only, a 
back slash ( \ ) indicates a period decrease only, a cross ( X ) indicates 
both increase and decrease of the period, and a filled circle ( • ) is used 
for systems for which have inadequate data for judgement. The magenta 
symbols are for the semidetached Algol-type binaries and the black ones 
are for detached Algol-type binaries. 



Figure 3. Plot of mass ratio(5) vs. secondary component's spectral type 
(SP2) for EW-type binaries listed in Table 3. The symbols are the same as 
in Figure 1 . 



0.1 and 0.8 were assigned to lower-precision observations (photo- 
graphic or visual ones) and high-precision observations (CCD or 
photoelectric ones), respectively. A weighted least-squares solution 
yields the following equation. 
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Min.I = 2427983.2814(±0.0188) + 4.''6296389(±0.0000035) x E 

-i-0.1977(±0.0096) sin[0.°0407 xE+ 163.°73(±0.°05)]. (2) 

The sinusoidal term in Eq. (2) suggests a periodic variation with a 
period of about 112.2 yr and an amplitude of about A = 0''.1977, 
which is more easily seen from the lower panel of Fig. 4, where the 
linear part of Eq. (2) was subtracted to the (O - C)i values. The 
good fit in Fig. 4 indicates no long-term steadily period increase or 
decrease. Therefore, we can exclude the presence of mass transfer, 
which is in accordance with the fact that WW Dra is a detached 
binary. The (O - C)2 values are shown in the fourth and eleventh 
column of Table 5. The residuals of the fit with Eq. (2) are dis- 
played in Fig. 5 and listed in the fifth and twelfth column of Table 
5. To detect possible regular trends in the residuals plotted in Fig. 
5, more high-precision times of light minimum are needed from 
future observations. 



Figure 2. Plot of mass ratio(^) vs. secondary component's spectral type 
(SP2) for EB-type binaries listed in Table 2. The symbols are the same as in 
Figure 1. 



filled circles to CCD or photoelectric ones. As shown in the upper 
panel of Fig. 4, the general (O - C)i trend can be described by a 
linear curve with superimposed a periodic fluctuation. Therefore, a 
sinusoidal term was added to a linear ephemeris to get a good fit 
to the (O - C)i curve (solid line in the upper panel of Fig. 4). To 
obtain a more accurate result, we focus the fit to only primary min- 
ima, though secondary minima also follow a similar trend. Weight 



5 DISCUSSIONS AND CONCLUSIONS 

In Section 4, we displayed the existence of a cyclical period change 
of WW Dra. This cyclical variation may be interpreted as due to ei- 
ther the magnetic activity of one or both components ( ApplegateJ 
1 1992) , or by the LTTE via the presenc e of a tertiary compan- 
ion. W ith the following equation given bv lRovithis-Livaniou et al. I 
( I2OOOI) , 



AP= V2[l -cos(2nP/P3)]xA, 



(3) 



where P^ is the period for the (O - C) oscillation, the rate of the 
period variation is found to be AP/ P = 5.29 xlO"'. In order to 
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Table 5. (O - C) data for WW Draco. 



JD.Hel. 

2400000+ 



(O-Oi (0-C)2 



JD.Hcl. 
2400000+ 



(O-Oi 



15205.599 

15501.809 

15774.848 

15955.582 

16390.638 

16418.486 

16603.702 

16969.44 

17131.500 

18867.661 

18955.495 

19191.586 

19890.777 

22881.486 

22895.458 

25247.247 

27284,275 

27284.32 

27307.396 

27321,296 

27335.197 

27344,447 

27534.261 

27543,530 

27557.409 

27557.418 

27645.3748 

27654.6307 

27691.6670 

27710.1918 

27881.4789 

27904.6277 

27918,5176 

27918.5197 

27932.4068 

27955.5536 

27983,3329 

28020.3700 

28057.4074 

28205.5510 

28219.4445 

28307.4054 

28404.6215 

33756.429 

33756.463 

34455.498 



-27607 
-2696 
-2637 
-2598 
-2504 
-2498 
-2458 
-2379 
-2344 
-1969 
-1950 
-1899 
-1748 
-1102 
-1099 



1247 
1247 
1398 



0.0619 
-0.0246 
-0.1339 
0.0444 
-0.0851 
-0.0148 
0.0158 
0.0128 
0.0357 
0.0844 
-0.0446 
-0.0649 
0.0515 
0.0177 
0.1008 
0.0364 
0.0260 
0.0709 
-0.0012 
0.0098 
0.0220 
0.0127 
0.0118 
0.0215 
0.0116 
0.0206 
0.0144 
0.0110 
0.0103 
0.0165 
0.0072 
0.0079 
0.0089 
0.011 
0.0092 
0.0078 
0.0093 
0.0093 
0.0097 
0.0050 
0.0096 
0.0075 
0.0013 
-0.0467 
-0.0127 
-0.0523 



0.27336 
0.18294 
0.07002 
0.24593 
0.11066 
0.18059 
0.20874 
0.2009 
0.22165 
0.24735 
0.11719 
0.09376 
0.2009 
0.12749 
0.2104 
0.11485 
0.07747 
0.12237 
0.04996 
0.06078 
0.07279 
0.06337 
0.05995 
0.06953 
0.05945 
0.06845 
0.06108 
0.05756 
0.05637 
0.06232 
0.05075 
0.05115 
0.05196 
0.05406 
0.05208 
0.05037 
0.05151 
0.05102 
0.05092 
0.04426 
0.04868 
0.04541 
0.03793 
-0.08096 
-0.04696 
-0.09582 



0.11874 
0.02288 
-0.09475 
0.07819 
-0.06367 
0.00585 
0.03144 
0.01897 
0.03786 
0.05087 
-0.07956 
-0.10356 
0.00341 
-0.04559 
0.03752 
-0.01305 
0.00209 
0.04699 
-0.02476 
-0.01355 
-0.00115 
-0.01031 
-0.00835 
0.00148 
-0.0082 
0.00079 
-0.00405 
-0.00731 
-0.00744 
-0.00095 
-0.00758 
-0.00652 
-0.0053 
-0.0032 
-0.00478 
-0.00581 
-0.00387 
-0.00329 
-0.0023 
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-0.00509 
0.03089 
0.06489 
0.03284 
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(1) 
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(1) 
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(1) 
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39029.556 

39205.511 

40844.335 

41043.396 

41154.547 

41168.406 

41682.288 

41682.298 

41763.327 

41830.474 

41904.533 

41904.536 

41918.503 

42617.379 

42904.6231 

42955.456 

43043.489 

43057.404 
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43161.535 
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43212.490 
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43307.414 
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43344.442 
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52576.270 
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54136.5301 

54210.598 
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54252.3029 



2386 
2424 
2778 
2821 
2845 
2848 
2959 
2959 
2976.5 
2991 
3007 
3007 
3010 
3161 
3223 
3234 
3253 
3256 
3259 
3278.5 
3284.5 
3289.5 
3293 
3310 
3315 
3318 
3329.5 
3415 
3496 
3556 
3556 
3556 
3648.5 
3737 
3929.5 
4244 
4655 
5109 
5234.5 
5277.5 
5312 
5515 
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5665 
5660 
5674 
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-0.0425 
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-0.1118 
-0.1191 
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-0.0766 
0.0011 
0.0426 
0.0428 
0.0430 
-0.1659 
-0.0224 
-0.0897 
0.1237 
0.1937 
0.1924 
0.1775 
0.2203 
0.3370 
0.3915 
0.4108 
0.4046 
0.4356 
0.4428 



-0.1756 
-0.14893 
-0.23664 
-0.25248 
-0.21415 
-0.24424 
-0.25834 
-0.24834 

-0.2225 1 
-0.23859 
-0.23559 
-0.15767 
-0.36553 
-0.16243 
-0.25621 
-0.18737 
-0.16146 
-0.21354 



-0.19113 
-0.15497 
-0.15937 
-0.16446 

-0.24381 
-0.17107 
-0.13325 
-0.13305 
-0.13285 



-0.0943 
-0.04955 
-0.07869 



0.05345 
0.0955 

0.10658 
0.0994 

0.13071 

0.13705 



0.01954 
0.04699 
-0.04031 
-0.05694 
-0.01914 
-0.04929 
-0.0666 
-0.0566 

-0.0319 
-0.04859 
-0.04559 
0.03219 
-0.18264 
0.01697 
-0.07745 
-0.00977 
0.01595 
-0.03631 



-0.01606 
0.01896 
0.01422 
0.00894 

-0.07735 
-0.01102 
0.0217 
0.0219 
0.0221 



-0.01512 
-0.02578 
-0.11832 



-0.01358 
0.00244 
-0.00261 
-0.01166 
0.02022 
0.02493 



)i used 
)i used 



)i used 



(1) 
(5) 
(6) 
(7) 
(7) 
(1) 
(1) 
(8) 
(1) 
(9) 
(9) 
(10) 
(11) 
(12) 
(1) 
(1) 
(1) 
(1) 
(1) 
(1) 
(1) 
(1) 
(1) 
(1) 
(1) 
(1) 
(1) 
(13) 
(14) 
(15) 
(14) 
(13) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 
(21) 
(22) 
(23) 
(24) 
(24) 
(25*) 
(26) 




o 
6 



o 
6 



0.4 












' 




r 


0.2 












. 


yi 


" ' 




? 


Oq^ 


o 


1 


o 
1 , 1 


^ 


1 


- 


0.0 
0.2 


<5 
1 


1 , 1 



-3000 


-2000 


-1000 





1000 


2000 


3000 4000 


5000 60 


0.4 






' 












0.2 
0.0 


■ o° 


1 o„ 


^o^>^. 


■^ 


^-^ 




«5*--^ 


^ 


0.2 


- 












-Tmm° 


- 


0.4 


- 


1 


, 




1 




o 


1 



























0.2 


















- 


CD 

w 

CD 


0.1 
0.0 


o 
o 

.°o« 


o 
o 





o 
o |B 


o 
oo 


^ oo o 


• 




0-: 


W 

CD 


-0.1 


o 
- o 




o 

















-0.2 


- 










o 
I.I. 






1 



-3000 -2000 -1000 



1000 2000 3000 4000 5000 6000 



-3000 -2000 -1000 



1000 2000 3000 4000 5000 6000 



Figure 4. (O - C) diagram of WW Dra. The upper panel : (O - C) i diagram 
of WW Dra computed witli Eq. (1). Tlie open circles refer to photographic 
or visual observations, filled circles refer to CCD or photoelectric ones. The 
soUd line refers to a combination of a Unear ephemeris and a cyclical period 
variation, and the dashed line to a new linear ephemeris. The lower panel : 
(O - C)2 curve of WW Dra as described by the sinusoidal term(solid line), 
after removing the linear term. The symbols are the same as in the upper 
panel. 



Figure 5. The residuals from fit witli Eq. (2). The symbols are the same as 
in the Fig. 4. 



reproduce this cyclic change, the required variation of the quadru- 
ple momentum AQ can be calculated with the following equation 
(Lanza & Rodono ,2002). 



AP 
~P 



Ma2 



(4) 



7 



where a is the separation between both components that can be 
determined with the Kepler's third law, 



Ml + M, = 0.0134 



p2' 



(5) 



where M, = 1.36Mo and Mj = 1.34Mo jAlbavraket aTl Il999l) . 
Combining Eq. (4) and Eq. (5), we derived AQi = 2.04 x 10'° g 
cm- and AQ2 = 2.01 xlO'" g cm^ for both components, respec- 
tively. Assuming conservation of the orbital a ngular momentum, 
the to tal AQ is on the order of 10^' - 10'^ g cm^ JLanza & Rodono I 
1 199^ ). which indicates the values of Agi and AQ2 for WW Dra 
are not typical ones for the close binaries, suggesting, that the 
mechanism of Applegate cannot interpret the cyclical period varia- 
tion of WW Dra. Moreover, the period of 1 12.2 yr for the {O - C) 
oscillation of WW Dra is too long in compariso n with magnetic 
cycle s in solar type sing le stars and close binaries JMaceroni et al. I 
Il99d : iBianchini 1119901) . Therefore, the cyclical period change is 
more plausibly interpreted as due to the presence of a third body. 

The good sinusoidal fit shown in Fig. 4 suggests that the orbit 
of the third body is approximately circular, which is different from 
the result derived bv lAlbavrak et al. I ( 119991) . Using ai2sin(' = Axc, 
where i' is the inclination of the orbit of the third component and c 
is the speed of light, ai2sini' is computed to be 34.25(±1.66) AU. 
Then combining the following well-known equation. 



4;r^ 



with 



/(m) = 



(M3 sin i'jf 
(Ml +M2 + Mif 



(6) 



(7) 



the mass function of the third body is computed to be f(m{) = 
3.19(±0.47)Mo. In the formula. Mi, M2, and M3 are the masses 
of the eclipsing pair and the third companion, respectively, G 
is the gravitational constant. According t o the same parameter s 
(Ml = 1.36Mo and M2 = 1.34Mq) used bv lAlbavrak et al. I fl999h . 
the lowest mass of the third body is calculated to be M3 = 6.43Mq, 
and the third body is orbiting the binary at a distance shorter than 
14.4 AU. When the third bod y is coplanar to the eclipsing binary 
: f = ; = 81. °4 (according to Albavrak e t al. ( 1999)), its mass is 
M3 = 6.57Mo. Using the formula given bV iMaver, (il99Qi) , 



Kr 



2k ai2sini3 



^3 VT 



(8) 



where Krv, P3, fli2 are in kilometer per second, years and AU, re- 
spectively, and considering the simplest situation of (3 = 90°, the 
semi-amplitude of the system velocity accompanied by the light- 
time eifect is approximately calculated to be 9.09 km s~ ' , which is a 
little less than the value determined by Albavrak et al. I J1999). Ac- 
cording to Allen's tables JDrilling & Landolt 120000 , the third com- 
panion is estimated to be a ~ B4 star. Therefore, it could be dis- 
covered by spectroscopic observation. However, no spectral lines 
of the third body were discovered up to now. It may be explained in 
two possible ways: (1) the star was observed in the past in a spec- 
tral range where the third body has no lines, or lines were present 
but the poor resolution of available spectra did not allow to detect 
them. Actually, it is difficult to find sufficient spectral lines to deter- 
mine radial velocity of B stars because their rapid rotational veloc- 
ity makes them too broad and weak to be accurately measured, or 
(2) the third body is a candidate black hole and it may play an im- 
portant role in the evolution of this system. The situation resembles 



that of V Pup J Oian et al. l2008a) . More observations are needed to 
check this hypothesis in the future. All these make WW Dra a very 
interesting system to study. 
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